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Abstract 

A cell dynamics method for domain separation of diblock copolymers (DBCPs) interacting with 
nanoparticles (NPs) whose diffusion coefficients depend on chain configuration is proposed for self- 
assembly of DBCP/NP composites. Increasing NP concentration slows down domain separation, 
but matching NP diffusion lengths and lamellar size of DBCPs reduces this effect. The model 
also explains features of different nanocomposites, such as morphological transitions induced by 
NPs, the coexistence of lamellar and hexagonal patterns in a single sample and peaked NP density 
profiles across the parallel domains. 
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I. INTRODUCTION 



Mixtures of diblock copo 



interest in recent years 



mers (DBCPs) and nanoparticles (NPs) have attracted much 
l|, |2j. The microphase separation of the DBCPs is expected to 
template NP arrangement, which may improve the physical properties of the composite or 
facilitate the production of other nanostructures of technological interest. However, the 
interaction of DBCPs and NPs may lead to a nontrivial morphology even with low NP 
loadings , as anticipated by models of mixtures in thermodynamic equilibrium 
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This motivated intense theoretical work 



10] , usually focusing on equilibrium properties. 



However, the most highly ordered patterns are obtained in far from equilibrium conditions 
rj ] and typical self-assemble times are of several hours or days. Thus, understanding the 
pathways to form the desired structures is essential to improve production methods. 

This paper introduces a model for self-assembly of a mixture of DBCPs and ex-situ formed 
NPs that describes features of various real composites. A cell dynamics method (CDM) 
I2I . represents the main physico-chemical mechanisms of domain separation of DBCPs 
and provides a realistic nanopattern morphology. The CDM variables, which represent the 
local chain configurations, interact with NPs whose diffusion coefficients depend on those 
configurations. Our CDM/NP model resembles those of Refs. Q, Ijj], but it is a significant 
extension of them because it accounts for a wider range of polymer-particle interactions. 
Compared to other recent models, the advantage of the CDM/NP approach is to address 
non-equilibrium and steady state features simultaneously. 

Among the non-equilibrium features represented by our model, emphasis is given to the 
role of particle mobility on the slow down of microphase separation and the effects of NP con- 
centration. High NP mobility is shown to destroy phase separation, while very low mobility 
delays the alignment with an interacting surface by freezing internal domain configurations. 



Connections with experimental work 
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17| are discussed and acceleration of global or- 



dering is suggested to occur when the NP diffusion lengths at the pure DBCP ordering 
time matches the lamellar size. The model also explains equilibrium features observed in 
experiments, such as the peaked or broad distributions of NP positions across DBCP lamel- 
lae [l9 1, the morphological transitions induced by NP loading 
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and the 



coexistence of different domain patterns in a single sample interacting with a surface 



fll 



The rest of this work is organized as follows. In Sec. II, we present the model and 
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discuss the simulation procedure. In Sec. Ill, we discuss the non-equilibrium evolution of 
the DBCP/NP composites, with a focus on the case in which the mixture interacts with an 
aligning surface. In Sec. IV, we show how the increase of NP concentration may change the 
morphology of the composites, in some cases producing samples where different structures 
coexist. In Sec. V, we show the relation of NP diffusion coefficients and the distribution of 
NPs across the ordered domains. In Sec. VI we summarize our results and conclusions. 



II. MODEL AND SIMULATION PROCEDURE 



In the CDM, a field variable ip (r, t) at cell f represents the state of a nanoscopic region 
of a DBCP sample at time t. The cell size is the length unit of the model and this is also 
the particle size, as discussed below, thus it is expected to correspond to lengths between 
2nm and lOnm in most applications (may be slightly larger in some cases). We will work 
with lamellar sizes which are typically ten times larger than the cell sizes, which in most 
cases correspond to high molecular weight polymers. 

In a unit time interval, i/j evolves as 

V (r, t+l) = ((1> (f, *)» + T (f, t)-((r (f, *))> - Bij (f, t) + £ (f, t) , (1) 

where T contains short-range interaction terms, —Bip (r, t) accounts for long-range inter- 
actions that restrict domain coarsening to a nanoscopic lengthscale, and £ is a thermal 
(conservative) noise. AX = ((X)) — X is the discrete isotropized laplacian of X. Here we 



follow the prescription of Ref. 24j for simplifying the short-range terms. In the case of 
symmetric DBCP, we choose 

T (r, t) = fty (f, t)] + Vhip (f, t) - bA 2 ^j (f, t) , (2) 

with 

/(^) = (1 + t)V-t^ 3 . (3) 

The first contribution to T represents attractive A-A and B-B interactions that lead to seg- 
regation of blocks A and B of the DBCPs, since / is a mapping with a flow to two symmetric 
stable points (it may be derived from a ip 4 approximation of a free-energy functional). The 
second term in Eq. [2] is related to local diffusive dynamics and the third one accounts for 
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;he curvature energy of the AB interface. For a review on the model and results, see Ref. 

CDM is a simplification of the physico-chemical processes of coarsening systems: the 
interaction terms have a simple form just to push ip to the stable points, the same diffusion 
coefficient is assumed for all chain configurations etc. Thus, it represents universal scaling 
properties of DBCP domain separation, particularly those morphological properties mea- 
sured on lengthscales near and above the lamellar size, but CDM is not expected to predict 
reasonable values of microscopic quantities such as local chain densities. 

Most works assume that if} is a difference of volume fractions of A and B ((f) a — 4>b) in 
a cell. However, this assumption makes typical CDM results inconsistent with properties of 



strongly segregated DBCPs (SSDBCPs). For instance, Fig, la s 
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rows the -0 profile obtained 
241 ] . It is clear that is 



24fl). Instead, SSDBCPs have 

1 in 



from CDM using typical parameters of the literature 
well below 1 in the middle of some domains (see also Ref. 

large regions with 100% A or 100% B separated by narrow interfaces, thus \<$a — 4>b 
the middle of those domains. One idea to overcome this problem is to work with parameter 
ranges corresponding to much lower temperatures, e. g. large values of r (see Refs. 24], [25]). 
However, this produces frozen configurations with very small domains (nearly the cell size) 
for extremely long times. 

On the other hand, the pattern morphology of real SSDBCPs is well represented by 
separating domains of positive and negative in CDM lattices shown in the literature - e 
g. compare images in Refs. []]• 
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24j . This feature 



|26( and CDM patterns in Refs. l3, 
and the universality character of CDM (i. e. the fact that CDM always includes the basic 
physico-chemical mechanisms of the coarsening process of SSDBCPs) lead us to propose a 
different interpretation of the field variable: (i) the sign of ip represents the type of chain 
in a cell, positive for A and negative for B; (ii) the absolute value of tj) is related to the 
local density at the cell and characterizes a certain position across the domain (ip may be a 
non-monotonic function of the local density, to be defined for a specific application); (iii) a 
range of small values of (as defined below) characterizes interfacial and/or mixed regions. 

The NPs are assumed to occupy approximately the whole volume of a cell. The state 
ijj of the CDM cell containing a NP represents neighboring chains interacting with it, i. e. 
the NP environment. The NP surface chemistry and this environment are responsible for 
the creation of energy barriers for the NP to move III, |27| . The NP diffusion is assumed 
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to take place in nearly equilibrium conditions with the local environment: the higher the 
barrier, the lower the NP mobility, which indicates a more stable configuration for NP 
localization. Thus, a NP at position r executes random walks among neighboring cells with 
a diffusion coefficient D [ip (r,t)], i. e. the rate with which the NP moves to neighboring 
cells is determined only by the chain configuration at the cell it moved from (and, of course, 
by the NP properties) . Target cells for the NP steps are randomly chosen among the nearest 
neighbors. Particle-particle interaction is restricted to the excluded volume condition. 

However, the presence of a NP in a cell also affects the states of neighboring chains. 
It is expected that those states evolve not only according to the basic (symmetric) CDM 
mechanisms, but also that neighboring chains move towards a certain configuration ipp that 
best wets the NP. For this reason, at each time step the state of the NP cell is replaced by 

^p(t + l) = (W)+g(^p-(W)), (4) 

where the first term represents the current neighborhood of the NP and the second one plus 
the CDM rules represent biased chain diffusion towards the state ipp. In this context, q 
measures the bias intensity. Following this reasoning, the states of NN and next NN (NNN) 
cells change as 

TpNN -> ?pNN + q/6 {tpP ~ IpNN) (5) 

and 

IpNNN -> tpNNN + g/12 (lp P ~ 4>NNN) , (6) 

respectively (conservation is ensured by suitable changes in their neighborhood). 

Although the parameters D, q and ipp were interpreted separately, they are related to 
the same microscopic NP-DBCP interaction, which is determined by DBCP properties and 
surface chemistry of NPs. For a given application, those parameters cannot be viewed as 
independent variables. However, the diversity of possible interactions in different systems 
and the difficulty to predict the real interaction terms suggest to explore them as indepen- 
dent for qualitative applications (in cases of solvent annealing, those parameters may also 
represent interactions with the solvent). For comparison we recall that previous CDM/NP 
models assumed that NPs were permanently embedded in a ib = 1 shell and had constant 

nn 

diffusion coefficient [14], H5(, thus showing the extension of the present approach. 

Our simulations were done in square lattices, so that relatively large systems could be 
simulated up to long times with a reasonable computational effort (this is necessary due to 



the large number of parameters to control). A disordered configuration at t = is represented 
by random distributions of ip in the range [—0.025,0.025]. In some cases we used periodic 
boundaries, but a large fraction of our simulations use a fixed boundary condition (FBC), 
where the row y = always has ip — 1, in order to represent a material preferentially wet 
by block A. 

A suitable set of CDM parameters (without NPs) were chosen to give the typical ip 
profile shown in Fig. lb after a domain separation time of order tsep ~ 10 5 . This large 
value of t S ep is interesting because it allows the study of the interplay between time scales 
of different processes. In Fig. lb, the cells closest to the A-B interface have —0.1 <ip< 0.1, 
thus this interval is associated to interfacial and/or mixed regions. The middle of A and B 
domains are characterized by |^| ~ 0.3 — 0.4 and the lamellar width is A ~ 10. The domain 
separation pattern obtained under these conditions, shown in Fig. lc, is similar to those 
presented in the literature, typically obtained with the parameters of Fig. la. 

Five simple models for D (ip) are used to investigate effects of NP mobility. 

Letting Da, Db, and Dj stand for A, B, and interfacial regions, respectively, the first 
three models are illustrated in Fig. 2a and defined as: DM1, with Db = \, Di = |, Da = §; 
DM2, with D B = f, D 1 = i D A = £q\ DM3, with D B = ^ D r = ^ D A = ^ 
(these values are given in time and length units of CDM defined above). DM1 represents 
large mobility in both domains, while DM2 and DM3 have Db/Da ~ 10 2 — 10 3 [as a rough 
guide, these are the ratios for Ag NPs in the domains of poly(styrene-b-methylmethacrylate) 
(PS-PMMA) of Ref. [ll|]. This form of D (ip) is useful to test the effect of different particle 
mobilities in different domains independently of possible intra-domain variations of D(ijj). 
With these models, any value ipp > 0.1 is reasonable for representing the effects of NPs on 
the chains because NP mobility is uniformly low in the whole A domain, where ip > 0.1. 

Fig. 2b shows the other two diffusion models, where D (ijj) is significantly reduced in 
certain chain environments, simulating deep energy valleys for NP diffusion at those po- 
sitions. For consistently describing particle-polymer interaction, the preference of the NP 
for a certain state [minimum of D (ijj)] must be accompanied by a flux of the neighboring 
chains towards that state. Thus, models DM4 and DM5 are respectively accompanied by 
ipp = 0.35 (middle of the A domains - see Fig. lb) and ipp = 0.1 (A side of the interfaces). 
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III. NON-EQUILIBRIUM SELF-ASSEMBLY OF THE COMPOSITE 



First we consider the effects of NP mobility on the global alignment of the composites 
with the FBC. 

In Fig. 3a, we show a FBC configuration without NPs obtained at t = 10 6 , which has 
ip profiles similar to Fig. lb except very close to y — 0. In Figs. 3b-d we show FBC 
configurations at the same time with 0.4% of NPs and diffusion models DM1, DM2 and 
DM3, respectively. The effect of NPs on the DBCPs is strong in this case, i. e. q is large. 
It is clearly observed that these small NP concentrations can slow down domain separation 
or eventually suppress the ordering. 

The quality of the alignment with the y = border is characterized by the correlation 
function along the y direction, 

C y (y, t) = (V (x', y', t) ^ (x\ y' + y, t))/([t/, (x\ y', t)f>, (7) 

with the average taken over x', y' and various realizations. The linear fits of the maxima of 
log \C y (y,t)\ x y plots give the corresponding correlation lengths L y . Fig 4a shows the time 
evolution of L y in the same conditions of Fig. 3a-d. Increasing L y indicates better alignment 
of the DBCP lamelles with the interacting border, i. e. a faster global organization of the 
composite. From Fig. 4a, it is clear that the best alignment with NPs is obtained with 
diffusion model DM2, followed by DM3 and DM1. 

In order to understand these features, an useful quantity is the diffusion length of NPs 
inside block A at a certain time t, which is given by 

l D = (D A t) 1/2 . (8) 

It measures the typical distance swept by a NP after that time if it is confined to the A 
domain, where it has the lowest mobility. 

For DM1, L y seems to saturate at long times (Fig. 4a), indicating that only a finite region 
near the bottom surface is aligned, while the rest of the sample is disordered. Simulation 
without the aligning surface provides a completely disordered lattice, i. e. there is no phase 
separation even at long times, ifjp = 0.35 and q — 0.1 were used in those simulations, so 
that the A chains are rapidly pulled to the NP neighborhood, but the fast movement of 
the NPs do not allow the chains to reorganize around them. Indeed, at the separation time 
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tsep of the pure system, we have Id ~ 112, which is much larger than the lamellar size 
A and explains the disordering effect of NPs. These NPs are equivalent to a thermal noise 
that brings the system above an order-disorder transition (ODT). It parallels the finding of 



Ref. 



161 ] that the ODT temperature of poly(styrene-b-isoprene) (PS-PI) decreases as the 



concentration of mixed silica particles increases. 

On the other hand, very small mobility (DM3) does not disturb phase separation in the 
middle of the lattice, as shown in Fig. 3d. Indeed, fixed NPs may improve local organization, 
as discussed in Ref. 28J. However, they slow down the alignment with the border y = 
because they may freeze a local order which does not match the order near the bottom 
surface. The slowly increasing global ordering is confirmed by the slow increase in L y of Fig. 
4a. The NP diffusion length inside the A domain at t$ep is Id ~ 2, while Id ~ 7 at t = 10 6 
(Fig. 2d). These values of Id are both smaller than the lamellar width A ~ 10, justifying 
the interpretation of these slow particles as pinning centers of a local domain configuration. 
The internal domain separation much faster than the global ordering of the sample was 
experimentally observed in composites of poly(styrene-b-2-vinylpyridine) (PS-P2VP) and 
CdSe particles [17]. 

The fastest global ordering (largest L y ) is obtained with DM2, where NPs in block A have 
Id ~ 22 at t$ep- This value of Id is near the domain period 2A, thus NPs can solidarily move 
with the neighboring chains to the most favorable positions without generating undesired 
noise. Other values of q and ipp lead to the same result. Thus, matching NP diffusion length 
and lamellar period at the pure DBCP ordering time t$ep may reduce the organization time 
To of the composite with the same morphology. In general, increasing temperature increases 
D and tsep, thus it also increases Id (but may decrease r c by bringing it closer to r SEP ). 
This may be crucial for technological applications where the experimental separation times 
upon annealing are several hours or days. 

As expected, the effect of increasing NP concentration is to increase the time for global 
alignment, or reduce L y at a fixed time. This is shown in Fig. 4b for diffusion model DM3 
(very small NP mobility). It parallels the increase of the characteristic time for lamellar 



23]. 



alignment of PS-PMMA films with Ag particles reported in Ref. 

In Figs. 4a and 4b, we note that L y scales diffusively with time (i. e. as t 1//2 ) in the pure 
DBCP, but the presence of NPs seems to lead to slower coarsening. It is not clear for us if 
there is a change in the universality class of this non-equilibrium system due to the NPs or 
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if this is a consequence of corrections to the main scaling form. This may be an interesting 
point for future theoretical investigation, but certainly much more accurate data in larger 
system sizes would be necessary. 



IV. MORPHOLOGICAL CHANGES 



The increase of NP concentration also has drastic effects on the composite patterns at 
long times. The steady states of the CDM/NP model are able to represent experimental 
features of some composites in thermodynamic equilibrium. 

First we analyze the effect of NP concentration with no aligning border. Configurations 
at t = 10 6 are shown in Figs. 5a-c, respectively without NP and with 0.8% and 3.2% of 
NPs. The pure DBCP has highly parallel lamellae (Fig. 5a) and this structure has only 
small deformation with low NP concentration (Fig. 5b). However, it changes to a hexagonal 
structure with 3.2% of NPs (Fig. 5c), even with weak effects of the NPs on the chains 
(q = 0.02). This is not a lattice effect because the CDM used square grids. This is a 
morphological transition caused by the stretching of the A domains that wet moving NPs, 
and can be achieved with smaller concentrations for larger values of q. 

Such transitions were experimentally observed in Ref. [20j with PS-covered Au particles 



in PS-P2VP. In Ref. 



211 ] . metallic and magnetic NPs changed the morphology of PS-PI 



from cylindrical to spherical domains; the two-dimensional projections of those patterns 
correspond to the lamellar-hexagonal transition illustrated in Figs. 5a-c (see Fig. 2 of 



Ref. 



2l|). Refs. [20] and 



22j show the opposite transition, from hexagonal to lamellar 



morphology, but these are also cases where NPs stretched one of the domains. Finally, the 
hexagonal morphology of PS-P2VP films was preserved with high loadings of CdSe particles 
in Ref. 29J], which suggests that suitable NP coatings were able to reduce NP effects on the 
chains. 

Even before these experiments were performed, NP-induced morphological transitions 



were theoretical 
Refs. a a 20 



predicted using different thermodynamic equilibrium models, e. g. in 
|. The stretching of A domains by NPs was already shown in previous 
CDM/NP models fl^ . \\\ . but no morphological transition was reported there. A recent 
work with dissipative particle dynamics (DPD) 31[ only showed deformation of the lamellar 
structure of DBCP with NPs. Thus, the novelty of the present approach is to show how 
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those transitions arise from a non-equilibrium cooperative dynamics. 

The increase of NP loading may also lead to the coexistence of different patterns in a 
single sample in contact with an aligning surface (FBC). This is shown in Figs. 6a and 6b, 
where the samples have parallel lamellae near the surface and nearly hexagonal order far 
from it. The increase of NP concentration and the preferential location near the interfaces 
(i. e. ipp ~ 0.1) favor the formation of the hexagonal structure, since the NPs have larger 
diffusion lengths near the interfaces and more easily stretch them. The concentration to 
obtain the coexisting patterns may vary, but a minimum cooperative effort of the set of NPs 
(high q or high concentration) is required to produce it. Similar structures were not shown 
in previous theoretical work, but it is interesting to recall that Ref. [32| showed evidence 
that non-selective NPs (which tend to go to the interfaces) favor morphological transitions 
from a lamellar phase, as illustrated here. 

Different patterns in a single sample were already observed experimentally. Thick films 
of PS-P2VP containing PS-covered Au particles 5] showed a lamellar structure near the 
air-polymer interface and a cylindrical structure (cross section with hexagonal packing) 
near the substrate. The lamellar region propagated to longer distances compared to Figs. 
6a-b, possibly affected by solvent evaporation and reduction of local NP concentration - 
similar effects of solvent concentration were reported in Ref. 33j. Moreover, in Ref. 231 ] . 
different patterns in a single sample were clearly observed: large Ag particles froze lamellar 
configurations of PS-PMMA perpendicular to the aligning (but relatively distant) substrate. 
The particular NP location in that case was related to the interaction with the air interface, 
but the overall effect of NPs with low mobility is similar to that in Fig. 3d. 



V. NP DISTRIBUTIONS ACROSS ALIGNED DOMAINS 

The NP position distributions across the parallel lamellae (FBC) can be also related to 
NP diffusion. 

If the NP diffusion coefficient is not very different in A and B domains (model DM1), 
then large fractions of NPs are found in the B domains. For instance, for low values of q (e. 
g. q = 0.02), a parallel lamellar organization is obtained with more than 10% of NPs in the 
B domains. 

However, a significant localization of the NP inside the A domains is obtained for models 
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DM2 and DM3, where the mobility in the A domains is much smaller than in the B domains. 
In Fig. 7a, we show long time probability densities P (y) of finding a NP at position y 
(measured relatively to the center of a lamellae) for DM2 with different values of ipp. In 
both cases the NPs are distributed almost uniformly inside the A domains, with sharp 
concentration decrease at the interfaces. Usually less than 1% of the NPs are in the B 
domain and less than 15% in the interfacial region. 

Similar shapes of P (y) are found for a broad range of parameters q and ipp and models 
DM2 and DM3. Thus, if the overall effect of surface chemistry and polymer properties 
is to reduce NP mobility across the whole A domain, then the NP distribution is broad, 
independently of the chain configuration that tends to move to the NP neighborhood. In 
these cases, the adjustments in chain configuration around the NPs (bias towards the state 
ipp) are balanced out by thermal effects, while slow NP movement avoids undesirable noise. 

On the other hand, highly peaked P (y) are obtained with DM4 and DM5, as shown in 
Fig. 7b. These are the cases where the NPs have severely reduced mobility in a certain 
position of the A domains [minimum of D (-0) - Fig. 2b]. This position determines the 
preferential localization of NPs: the center of the A domains for DM4 and the interfaces for 
DM5. 

These results suggest that the localization of NPs strongly depends on the surface chem- 
istry being able to provide a favorable interaction with a particular chain configuration, 
which may be the middle of the A domains or the interface in the above examples. This 
favorable interaction corresponds to some energy valley at that position, which reduces NP 
diffusion coefficient. Of course additional conditions are important for the global organiza- 
tion, such as reduced NP effects on the chains and avoidance of pinning effects due to very 
low NP mobility, as discussed in Sec. III. 

The above NP distributions resemble those of Kim et al 



18 



19| for Au particles in thick 



films of PS-6-P2VP: preferential localization of the NPs in the center of the PS lamellae or 
at the interface with P2VP, depending on the density of coatings (PS or P2VP) of the NP 
surface. Recently, the same group [34| showed NP position distributions similar to those in 
Figs. 7a-b using various coatings. This suggests to search for an interpretation based on 
position-dependent diffusion energy barriers during the annealing process. 
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VI. SUMMARY AND CONCLUSION 



We proposed a cell dynamics model for DBCP domain separation interacting with ran- 
domly moving NPs. Increasing NP concentration slows down separation and global align- 
ment, but this effect is reduced by matching the lamellar size and NP diffusion lengths. The 
model also represents various features of real composites, such as morphological transitions 
induced by NPs, the coexistence of different patterns in a single sample and various NP 
density profiles across the lamellae. 



35, 
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381 ] and mor- 



Previous works alsopredicted localized NP distributions 
phological transitions jg, 0, [^(J in thermodynamic equilibrium, but an advantage of the 
CDM/NP approach is to explore the non-equilibrium routes leading to those and other 
properties. An additional advantage is the low computational cost; for instance, one config- 
uration of a 256 x 256 lattice at t — 10 6 is generated in less than two hours in a desktop. 
This is important if one aims at testing hypothesis theoretically before performing series of 
experiments. Quantitative application is possible by associating the time and length scales 
involved in the domain coarsening ( e. g. t$ep and A) to experimental values and choos- 
ing sensible values of the other parameters. Certainly the model can also be extended, for 
instance to include in-situ NP formation or NP interaction. 
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FIG. 1: (a) Long time ip profiles in directions perpendicular to the parallel DBCP domains without 
NPs and traditional parameters of the literature: r = 0.3, D = 0.5, B = 0.005, b = C = 0. 
(b) Profiles with typical parameters of this work: r = 0.035, D = 0.5, B = 0.0005, b = 0.35, 
C = 0.02; the region between dashed lines is the interfacial and/or mixed region, (c) Typical 
lattice configuration at tsep = 10 5 obtained with the CDM parameters of this work. A (B) 
domains are light gray (white) and interfacial regions are medium gray. 
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FIG. 2: Diffusion models: (a) DM1 (dashed lines), DM2 (solid lines), and DM3 (dotted lines); (b) 
DM4 (solid line) and DM5 (dashed line). 
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(a) (b) 




FIG. 3: Typical lattice configurations at t = 10 6 with the FBC: a) without particles; b), c), d) 
0.4% of NPs with DM1, DM2 and DM3, respectively, q = 0.1, and ip P = 0.3. Surface with fixed 
■0 = 1 is dark gray, A (B) domains are light gray (white), interfacial regions are medium gray, and 
NPs are black. 



17 




FIG. 4: Time evolution of correlation length with the FBC: (a) same conditions of Fig. 2, with no 
particles (squares), DM1 (triangles), DM2 (crosses), and DM3 (asterisks); (b) diffusion model DM3 
with the same q and ipp of Fig. 2, and NP concentrations 0.4% (asterisks), 0.8% (filled triangles), 
and 1.6% (filled squares). Error bars are near the size of the points. 
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(a) (b) 




FIG. 5: Typical lattice configurations at t = 10 6 without (a), with 0.8% (b) and with 3.2% (c) of 
DM2 NPs, no aligning boundary, b = 0.05, q = 0.02, and ipp = 0.3. Color scheme is the same of 
Fig. 2. 
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FIG. 6: Typical lattice configurations with FBC at t = 10 7 and: (a) 0.8% of DM2 NPs, b = 0.35, 
q = 0.1, and = 0.1; (d) 3.2% of DM5 NPs, b = 0.35, q = 0.02, and ^ P = 0.1. Color scheme is 
the same of Fig. 2. 
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FIG. 7: NP position distributions across the parallel lamellae for: (a) DM2 with tpp = 0.1 (solid 
line) and i> P = 0.35 (dashed line); (b) DM4 with i\) P = 0.35 (solid line) and DM5 with ip P = 0.1 
(dashed line). Horizontal bars indicate typical interface regions. 
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